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PREFACE 


A  similar  report  was  issued  30  December  1957  at  Space 
Technology  Laboratories,  Inc.  ,  "Optimum  Propellant 
Loading  for  Titan  and  Thor  Missiles"  (Confidential), 
by  David  W.  Whitcombe.  In  the  present  unclassified 
edition,  Thcr  and  Titan  performance  data  are  excluded. 
Sections  have  been  added  describing  the  calculation  of 
mean  outage  and  the  operation  of  a  propellant  utili  ~tion 
system. 


ABSTRACT 


This  report  derives  optimum  fuel  or  mixture  ratio  biasing  techniques 
commonly  employed  in  current  ballistic  missile  programs.  The  missile  stages 
can  use  conventional  loading  ,  g.  ,  Titan  and  Thor)  or  they  can  use  propellant 
utilization  (PU)  systems  (c.  g.  ,  Atlas  and  Centaur).  A  description  of  the  mecha¬ 
nization  errors  leading  to  unburned  propellant  (outage)  is  given  for  both  cases. 
Formulae,  derived  in  the  report,  allow  the  calculation  of  mean  outage,  0N;  the 

2  2 

mean  square  outage,  <9  >;  the  outage  variance, o-^;  as  well  as  the  probability 
that  the  outage  is  less  than  some  fixed  value,  P(G  <  6^). 

The  analysis  in  this  report  applies  directly  to  conventionally  loaded  bal¬ 
listic  missiles.  It  is  assumed  that  the  missile  stage  is  loaaed  in  accordance 
with  a  loading  mixture  ratio,  r^.  That  is,  the  piupellant  is  loaded  such  that 
the  ratio  of  oxidizer  to  fuel  is  rL  in  the  nominal  missile.  A  mixture  ratio 
bias,  Ar,  is  calculated  that  will  minimize  0^  and  <ez>,  and  maximize 
P(0  <  0Q).  An  identification  is  obtained  that  extends  the  basic  loading  formulae 
to  stages  employing  PU  systems.  Additional  formulae  are  given  that  allow 
calculation  of  a  fuel  or  oxidizer  bias  equivalent  to  the  mixture  ratio  bias. 
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SECTION  i 


INTRODUCTION  AND  SUMMARY 

It  is  shown  in  this  report  that  the  total  propellant  per  stage  in  the  Titan 
and  Thor  missiles  is  best  loaded  with  the  aid  of  a  loading  mixture  ratio,  rT  . 

Li 

When  rT  is  obtained,  the  nominal  initial  masses  of  oxidizer  and  fuel,  M, 

L  LN 

and  Mp^,  respectively,  to  be  loaded  are  given  by  the  relationships: 


M 


M 


M 


N 


FN  "  TT 


LN 


MNrL 

TT~r7 


*(  3.2a) 


(3.  Zb) 


where  M^  is  the  total  initial  mass  of  propellant  in  the  nominal  missile.  The 
optimum  value  for  is  less  than  the  nominal  expected  burning  mixture  ratio, 
r^;  that  is, 


rL  =  rN  '  Ar 


(3.5) 


Formulae  for  the  coii»(Ja’..~t.:''*n  of  Ar  will  be  derived  in  this  report. 

Loading  the  Titan  or  Thor  in  accordance  with  r^  is  equi  -alent  to  loading 
the  missile  in  accordance  with  r^  and  then  adding  an  excess  of  fuel,  AMp,  or, 
equivalently  removing  a  portion  of  the  oxidizer,  AM^.  The  quantities  AMp 


''Equation  numbers  in  this  section  refer  directly  to  equations  derived  in 
Sections  2  through  8. 
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and  AM.  are  relatM  to  Ar  by  the  following  relationship: 


AM  .  -AM. 

F  Ar  L 

Tvi —  =  r~  ~  - 

FN  N  LN 


(3.  7) 


The  quantity  AMp  is  referred  to  as  the  "fuel  bias.  " 

An  attempt  is  made  to  load  the  Titan  and  Thor  to  the  desired  loading  mix¬ 
ture  ratio.  The  propellant  is  loaded  using  weight-measuring  load  cells  or  using 
volume -measuring  flow  devices.  In  either  case,  errors  in  quantity  of  both 
oxidizer  and  fuel  will  exist.  In  addition,  the  actual  burning  mixture  ratio  will 
differ  from  the  intended  burning  mixture  ratio.  As  a  result  of  these  errors,  an 
excess  of  either  oxidizer  or  fuel  will  remain  in  the  tanks  when  the  other  is 
exhausted.  The  residual  unburned  propellant  is  referred  to  as  outage.  Outage 
will  be  denoted  by  the  symbol  0.  and  will  be  always  be  expressed  as  a  fraction 
of  the  total  stage  propellant. 

A  propellant  utilization  {PU)  control  system  may  be  developed  that  will 
cause  the  propellant  to  be  burned  in  closed-loop  fashion.  Such  a  technique  is 
employed  in  the  Atlas  missile.  This  method  requires  measurement  during 
flight  of  both  the  remaining  oxidizer  and  the  remaining  fuel.  The  burning  mix¬ 
ture  ratio  is  then  adjusted  so  that  the  propellants  are  exhausted  simultaneously. 
Errors  in  the  measurement  of  remaining  o..idi/.er  and  fuel  result  in  outage  even 
when  a  PU  system  is  — '  The  outage  c«.i  !>■  mi'  imi?  d  by  introducing  a  known 

bias  in  the  fuel  measurement.  It  will  be  shown  in  lii,  '‘port  'hit  choosing  the 
fuel  bias  for  a  PU  system  is  analogous  to  choosing  the  loading  mixture  ratio  for 
the  Titan  and  Thor  missiles. 

The  following  statistical  properties  of  outage  are  derived  in  the  body  of 
this  report: 

.'9C>  -  me square  out.ug.. 

-  mean  outage 


P(9  <  9q)  =  probability  that  the  outage  is  less  than  some  fixed 
constant  (9 q) 


<Tg  =  outage  variance 


=  <e2> 


6 


2 

N 


<r 


9 


V 


Five  particular  methods  of  loading  a  missile  that  does  not  have  a  PU  sys¬ 
tem  are  discussed: 

1.  Standard  loading:  In  this  method  =  r^,,  the  nominal  expected 
burning  mixture  ratio.  This  procedure  is  discussed  in  Section  2. 

2.  Minimum  mean  square  outage  loading:  In  this  method  r^  =  r^  -  Ar, 
where  Ar  is  chosen  to  minimize  the  root  mean  square  (rms)  outage. 

The  method  for  determining  Ar  is  given  in  Section  3. 

3.  Maximum  P(9  <  6^)  loading:  In  this  method,  r^  =  -  Ar,  where 

Ar  is  chosen  to  maximize  P(0  <  6q).  This  method  is  discussed  in 
Section  4. 

4.  Minimum  mean  outage  loading:  In  this  method  rj_,  =  rjvj  '  Ar,  where 
Ar  is  chosen  to  minimize  the  mean  outage.  This  procedure,  dis¬ 
cussed  in  Section  5,  is  used  when  it  is  desired  to  maximize  the 

pel  formance  of  the  nominal  mis  >’le. 

5.  Maximum  P1'-  .  i  load:"?-  In  this  metl  disc  issed  in  Section  6, 

Ar  is  chosen  to  maximize  P(V  >  Vn),  wher/'  V  us  a  gone  perform- 
ance  parameter. 

An  analog  for  each  of  the  above  loading  procedures  exists  when  a  PU  sys¬ 
tem  is  employed.  The  only  change  is  that  a  fuel  measurement  bias, 

Af(  as  AMJ.  (rather  than  a  mixture  ratio  bias)  is  obtained  to  achieve  the 
r 

statfj  c c*sJ'c'v-.. 
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An  exarnpi':,  using  the  formulae  developed  in  this  report,  will  be  given. 

It  will  be  assumed  in  the  example  that  the  missile  does  not  have  a  PU  system. 
The  following  data  on  the  missile  will  be  used: 

Mfn  =  60,000  lb 

M.  ..  =  135,  000  lb 
LN 

Mn  =  195,000  lb 
rN  =  2‘  25 

The  following  data  on  propellant  loading  and  burning  wi1'  be  assumed: 

Initial  fuel  mass  loading  tolerance  (o-p)  =  0.  25% 

Initial  oxidizer  mass  loading  tolerance  (ir^)  =  0.  5% 

Burning  mixture  ratio  tolerance  (o^)  =  0.  5% 

In  Table  I,  the  quantities  e^,,v/<02>,  and  <i  g  are  give:;  as  fractions 
of  the  total  nominal  mass  of  propellant  (195,000  lb).  Note  that  when  the  mean 
outage  is  minimized,  the  rms  outage  is  near  minimum,  and  that  the  P(0  <  1%) 
is  near  maximum.  For  most  missions,  acceptable  results  .vill  be  obtained 
when  0j^j  is  minimized. 

A  similar  table  could  readily  bo  constructed  for  a  missile  stage  utilizing 
a  PU  system.  It  is  shown  in  Section  8  that  the  principal  change  in  the  analysis 
results  from  ;  <;dtfini.:g  «■  and  o-,  and  choosing  r 

Summary  oi  Equations 

The  more  important  relationships  and  definitions  derived  in  the  body  of 
the  report  will  be  summarized  in  the  following  paragraphs.  The  equations  will 
be  given  assuming  conventional  loading  procedures  as  employed  in  the  Titan 
'rv>.-\r  stages.  The  following  definitions  are  used: 

=  nominal  total  mass  of  propellant  loaded 


1-4 


s.  =  nominal  total  mass  of  fuel  loaded 


N 

1  t  rT 


(3.  2a) 


=  nominal  total  mass  of  oxidizer  loaded 

v*  r 
N  L 


(3.  2b) 


=  njminal  ;xpected  burning  mixture  ratio 

=  loading  mixture  ratio  =  r^  ■  Ar  (3.  5) 

=  initial  fuel  mass  leading  tolerances  ( I cr ) , 
expressed  as  a  fraction 

=  Initial  oxidizer  mass  loading  tolerance  ( lcr ) , 
expressed  as  a  fraction  of 

=  burning  mixture  ratio  tolerance  ( ! «r ) ,  expressed 

as  a  fraction  of  r.. 

h 

=  /4  +  +  %  (2-i5) 

=  (3.12) 


1  2 

1  ’?  ° 
<P'{o)  =  -L:  c 


=  me  :n  outage 


;0  >  -  mean  square  outage 


=  outage  variance  (about  the  mean) 
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aj  =  some  fixed  constant  value  of  outage 
P(8  <  0q)  =  probability  that  the  outage  is  less  than  0q 

/  2  ' 

The  quantities  8^,  sj  <8  >,  cr g ,  and  0q  are  all  measured  as  fractions  of  the 
total  mass  of  propellant. 

Mean  outage: 

The  expression  for  mean  outage  is  calculated  in  the  body  of  the  report  as 


=  P  0*  (n)  +  °  <t>  (a)  - 


a  rN 


N  '  1 
:N  +  1 


The  mean  outage  is  minimized  when 


,  ,  1  rN  ’  1 

P(t.)  =  77— n 


(5.  11) 


9m  =  P 0'  (°> 


(5.  12) 


■?  ? 
<e~>  =  p~ 


f  (1  +  o^)r„  r„,  -  1 


_(I  +  rN> 


7771 

N 


p-  i  n 


Table  III  may  be  used  when  the  propellants  are  loaded  to  minimize  the  mean 
cui  "<gi’  or  r-’i  ai  ■  aasticai  p?.-'.»neters  are  calculated  as  a  function  of  mixture 


ratio,  r 
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Mean  square  out?-;* e : 

The  mean  square  outage  is  obtained  as: 


2,  2  . 

2  1J  rN  *  2 

<6  >  = - - - 7-  (1  +  a) 

2(1  +  rN)^ 


r-,  -  1 


]  p2ja0'(n)  +  (1  +  a2)  0  (u)j  (3.  16,  5.  16) 


The  mean  square  outage  is  minimized  when 


<j>'  (a)  +  a0  (a) 


+  1 
-  1) 


a 


(3.  20) 


Then 


0 


N 


N 


<02>  =  o2 


N 


(1  +  rN)‘ 


l 

J 


The  outage  variance  may  be  determined  in  general  from  Equations  (3.  16) 
rnd  15.  5)  as 


»8  =  ^  -  9N 


(5.  13) 
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Probability  that  0  <  0Q: 


p(0  <  e0)  =  <p(*0)  +  ^(uj) 


(4.  11a) 


where 


..  .(1+rN>90. 

u  n  -  -  +  «a 

0  or 
H  N 


{4.  lib) 


I 


..  .  °  +  rN)60 

“l  -  — - a 


(4.  i  lc) 


The  quantity  P{0  <  0n)  is  maximized  when 


rN  '  1 

Q  =  TpTT"  90 

1  N 


(4.  14) 


In  this  special  case,  Equations  (4.  11)  become 


P(0  <  0Q)  =  20(uo) 


(4.  16) 


d  +  rN) 

"  2PrN  °o 


(4.  15) 
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Table  IV  may  be  used  when  the  propellants  are  loaded  to  maximize 
P(6  <  0q).  When  this  table  is  used,  it  is  assumed  that  0g/p  is  obtained  by 
substituting  Uq  =  3  into  Equation  (4.  15).  Several  statistical  parameters  of 
interest  are  calculated  in  this  table. 

Probability  that  6  =  fuel  or  oxidizer: 

The  probability  that  the  outage  will  be  all  fuel  or  all  oxidizer  ia  dependent 
on  a.  It  is  shown  in  the  text  that  these  probabilities  are  evaluated  as: 

P(0  =  fuel)  =  i  +  £(Q) 

P(0  =  oxidizer)  =  -  -  0(a)  (5.21) 


PU  System  Analysis 

The  following  definitions  are  used  in  the  analysis  of  PU  systems-. 

or c  =  rms  fuel  measurement  error  (l<r),  expressed  as  a  fraction 
ofMFN 


cr^  =  rms  oxidizer  measurement  error  ( 1  o' ) ,  expressed  as  a  fraction 

of  M... 

LN 


2 
*  i 


Af  =  fuel  measurement  bias  in  pounds  (to  b<.  -  ibtracted  from  each  PU 
system  fuel  measurement) 

All  outage  statistical  parameters  that  result  whet,  a  PU  system  is  employed 
may  be  obtained  using  the  relations  that  apply  for  conventional  loading  proce¬ 
dures.  It  is  only  necessary  to  make  the  following  identification: 


P  “•  q 


(8.  11a) 
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a  En 


Af 


(8.  lib) 


qMFN 

Note  that  the  PU  system  fuel  bias  Af  is  equivalent  to  the  fuel  loading  bias  dis¬ 
cussed  in  Equation  (3.  7a). 

The  quantities  <p(a)  and  <p'{a)  are  tabulated  in  Mathematical  Tables  from 
Handbook  of  Chemistry  and  Physics,  (Chemical  Rubber  Publishing  Co.  , 
Cleveland,  Ohio).  For  convenience,  an  abridged  tabulation  of  these  functions  is 
given  in  Table  II. 
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outage  0.2638  0.3344  0.9924  0.2055  0.0085  2.2415  226  609 


Table  II.  Tabulation  of  functions. 
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<P  (a) 

<p'  (a) 

<P'{  a)  +  a<p(a) 

0 

0 

0. 3989 

0. 39890 

0.  2 

0.  0793 

0.  3910 

0. 40686 

0.  4 

0.  1554 

0. 3683 

0.  43046 

0.  6 

0. 2258 

0. 3332 

0.  46868 

0.  8 

0. 2881 

0.  2897 

0. 52018 

1.  0 

0.  3413 

0. 2420 

0.  58330 

1.  2 

0.  3849 

0. 1942 

0.  65608 

1.  4 

0. 4192 

0.  1497 

0. 73658 
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0. 82322 
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C.  4°7/* 
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1. 40063 
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SECTION  2 


OUTAGE  (WITH  STANDARD  LOADING) 

In  this  section,  an  expression  for  the  mean  square  outage  for  a  single 
stage  is  derived;  the  calculation  involves  these  parameters: 

M  =  initial  mass  of  propellant 

M_  =  initial  mass  of  fuel 

Mr  =  initial  mass  of  oxidizer 

Lj 

r  =  mixture  ratio. 

Mass  flow  rates  are  indicated  as  Ivl,  Mp,  and  a  subscript  N  is  used  to 

indicate  nominal  values.  The  burning  mixture  ratio  is  defined  as 


M. 


Hence, 


w  _  M 
MF  "  T  +  r 


and  Mr  =  t 


•M 


i  H 


(2.  2) 


All  the  quantities  defined  above  are  assumed  to  be  constants  the  foil  awing 
discussion. 

The  missile  engine  will  operate  until  the  supply  of  either  fuel  or  oxidizer 
is  exhausted.  The  remaining  amount  of  fuel  or  oxidizer  is  called  outage.  Note 
that  outage  is  always  positive  and  is  always  all  fuel  or  jail  oxidizer.  The  outage 
csicu.  die'*  then  u  .pends  on  v.h.  mer  the  outage  is  all  fuel  or  all  oxidizer. 
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When  the  outage  is  all  fuel,  the  burning  time,  t^,  is  given  by 


t 


L 


(2.  3) 


The  mass  of  fuel  consumed  during  this  period  of  time  is  then 


•  ML 

AM„  =  M„t.  =  Mr<  - — 

f  f  l  r  ml 


(2.  4) 


The  fuel  outage  may  then  be  written  as 


8 


F 


-  M 


F 


(2.5) 


Similarly,  if  the  fuel  burns  out  first,  the  outage,  all  oxidizer,  can  be  calculated 
as 


8 


L 


\1 


! 


Mr 


M 
1  ! ' 


(2.6) 


The  above  expressions  may  be  written  as 


0F  =  MF  -  ±Ml  (2.  7a) 

0 l  =  Ml  -  rMp  (2.  7b) 
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Note  that  the  abc-a  two  events  are  mutually  exclusive  because  the  outage  is 
always  positive.  The  outage  is  calculated  using  Equation  (2.  7a)  or  Equation 
(2.  7b)  depending  on  which  gives  a  positive  result. 

There  can  be  no  outage  for  the  nominal  missile  in  the  standard  method  of 
propellant  loading!  that  is,  an  attempt  is  made  to  initially  load  the  missile  such 
that 


M 


M 


N 


FN  "  T~T 


and  M 


rNMN 


'N 


LN 


nr 


N 


(2.8) 


If  the  missile  were  actually  loaded  this  way  and  the  mixture  ratio  were  nominal, 
there  would  never  be  any  outage  in  a  nominal  missile.  However,  because  of 
ground  and  inflight  temperature  variations,  propellant  measurement  tolerances, 
and  burning  mixture  ratio  tolerances,  this  performance  can  never  be  achieved. 

The  linear  variation  in  outage  may  be  obtained  as  a  function  of  the  toler¬ 
ances  on  M„,  M.  ,  and  r.  This  result  is: 

X1  Lj 


50 


F 


6Mf 


+ 


6r 

~r 

1  TvT 


M 


LN 


(2.9a) 


or 


50l  =  6Ml  -  hr  Mfn  -  rN  6MF 


(2.9b) 


The  mean  square  outage  may  now  be  calculated  using  Equation  (2.  9a)  or  (2.  9b) 
which’ ve.  jjl.c-  a  'lusitivc  Assume  now  that  the  loading  is  accomplished 

according  to  Equation  (2.  8)  and  6MF>  6r,  and  6M.  are  uncorrelated  and 
normally  distributed  with  zero  mean.  Then  the  outage  is  equally  likely  to  be  all 
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fuel  or  all  oxidizer.  The  mean  square  outage,  <9  >,  written  as  a  fraction  of  tk' 
total  propellant,  may  then  be  evaluated  as 


2 

<9  > 


(2.  10) 


where  Xp  and  X^  are  mean  square  fuel  and  oxidizer  outages  expressed  as 
fractions  of  the  total  propellant.  These  quantities  may  be  calculated  from 
Equation  (2.  9)  by  squaring,  averaging,  and  dividing  by  the  square  of  the  nominal 
mass  of  propellant  as 


and 


v2  _(  (56F) 

Ap  -  55 


;> 


M 


N 


1  ,2.2,2. 

1  '°f  vl.  0  r' 


O  +  rN) 


(2.  ID 


M 


N 


2 

r 

N 


(1  +  rN) 


2  WF  +  +  <V 


(2.  12) 


2-4 


where 


(2.  13) 


and  <•  •  •>  indicates  an  ensemble  average.  It  is  assumed  in  the  above  equations 
that  the  missile  is  loaded  in  accordance  with  Equation  (2.  8). 

When  the  above  result  is  substituted  into  Equation  {2.  10),  the  following 
expression  for  mean  square  outage  is  obtained. 


r2  +  1 

<02>  =  '  I?-  -  .  p2 

(rN  +  l)2 


(2.  14) 


where 


<r„  =  rrns  fuel  tolerance  ()<r).  expressed  as  a  fraction  of  M™, 

I*  r  IN 

<t,  -  ms  oxidizer  tolerance  ( ltr) ,  expressed  as  a  fraction  of  M. 

i-j[N 

-  rms  burning  mixture  ratio  tolerance  (1<t),  expressed  as  a 

fraction  of  r^.  (2.  15) 
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SECTION  3 


MINIMUM  MEAN  SQUARE  OUTAGE  LOADING 

It  is  clear  from  Equation  (2.7)  that  the  expected  oxidizer  outage  is 
times  the  expected  fuel  outage  when  the  missile  is  loaded  m  the  standard  way 
[i.  e.  .  in  accordance  with  Equation  (2.  3)].  Thus,  an  attempt  should  be  made 
to  prevent  the  oxidizer  from  burning  out  first.  This  may  be  accomplished  by 
loading  more  fuel  and  less  oxidizer  than  standard,  keeping  the  same  total  mass 
of  propellant.  This  is  equivalent  to  loading  the  missile  on  the  basis  of  a 
loading  mixture  ratio,  r^,  where 


rL  <  rN 


Then  the  missile  would  be  loaded  such  that 


M. 


M 


N 


FN  1  +  rT 


M 


M..  rT 
N  L 


LN  I  +  r, 


(3.  1) 


(3.2a) 


(3.  2b) 


The  result  of  such  loading  is  .i  suiptus  ox  xucl  in  the  nominal  missile,  after  all 
the  oxidizer  is  burned.  Assume  now  that  the  missile  is  loaded  according  to 
Equation  (3.  2).  Then,  using  similar  methods  as  were  used  in  Section  2,  the 
linearized  fuel  and  oxidizer  outages  may  be  written  as 


60.,  =  6MC.  +  -~~  M 

J*  i* 

rN 


-  —  6M7  +  —  M_. 
LlvI  rN  L  rN  FN 


(3.3) 
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ana 


6  6r 


lN 


1  5Mt 


It.  mfn 

N 


6  M. 


LN 


At 

rN 


M 


FN  / 
/ 


(3.4) 


where 


Ar 


N 


(3.  5) 


To  the  first  approximation 


M 


LN 


N 


FN 


hence,  the  above  equations  may  be  written  as 


and 


6  0, 


M 


FN 


M 


N 


M 


N 


6  M_, 


6r 

r 


6M 


N 


M 


L  x  Ar 


N 


LN 


6eL 

'rNMFN  j 

UM 

6r 

6Ml 

+  Al 

mn 

IvIN  1 

,mfn 

rN 

mln 

rN 

(3.6a) 


(3.6b) 
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'*  6)  would  be  equally  valid  if  the  substitutions 


Note  that  Equation  !' 


Ar 


amf 


FN 


(3.  7a‘ 


or 


Ar 

rN 


-AM, 


M 


I.N 


(3.7b) 


were  made  in  Equations  (3.  6).  When  Equation  (3.  7a)  is  used,  the  fuel  bias, 

AM„,  rather  than  the  mixture  ratio  bias,  Ar,  would  be  obtained.  The  calcu- 
i 

lation  of  the  fuel  bias  is  convenient  when  it  is  desired  to  always  load  the  oxidizer 
tank  to  capacity.  The  fuel  tank  is  loaded  with  the  nominal  quantity  of  fuel  plus 
the  fuel  bias.  Equation  (3.  7b)  calculates  an  oxidizer  bias,  AM,  This  quantity 
would  be  useful  if  the  fuel  tank  were  always  loaded  to  capacity.  The  analysis 
used  in  this  report  will  be  presented  in  lerms  of  the  mixture  ratio  bias.  The 
fuel  bias,  or  the  oxidizer  bias  can  always  be  calculated,  ii  uosired,  using 
Equations  (3.7). 


Outage  may  now  be  calculated  using  Equations  (3.  6a)  or  (3.  6b),  whichever 
gives  a  positive  result  Equations  (3.6)  m«v  be  written  in  the  following  way, 


66  _  1 
MN  "  1  +  fN 


X,  X  >  0 


N 


>  +  r 


X,  X  >  0 


(3.8) 


N 
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where 


_  5Wf 

t  — 

6M. 

L  +  Ar 

mfn 

rN 

MLN  rN 

The  quantity  X  is  now  normally  distributed  with  mean  - —  .  Then  the  mean 
square  outage,  <  0^  >,  may  be  written  as  ^ 


62  > 


(1  +  rN) 


x  F(x)  dx  + 


2  r  0 

r..  i  _ 

- j  I  x  F(x)  dx 

( 1  +  fjq)  -  oo 


(3.  9) 


i 


where 


F(x) 


i 

—  ■  -  e 
\!  Zr,  p 


Ar\  2 


(3.  10) 


and,  as  defined  in  Equation  (2.  15), 


2 

P 


•> 


IT 


2 

r 


Equation  (3.  10)  may  now  be  substituted  into  Equation  (3.  9).  After  a  change  in 
variables,  the  result  is 
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2 

<  e~> 


0  +  rN) 


0'  (t)  dt 


+ 


2 

rN 


(1  +  rjq' 


2  <*>' (t)  dt 


where 


*'(t) 


2 


t 


2 


a 


Ar 

PrN 


Expanding  the  parentheses  in  Equation  (3.  11)  and  simplifying  gives 


<e2> 


I  y 

I  t  *’(t)dt 
i 

J- a 


r°°  . 

2n  /  t  ys'  (t)  dt  +  a-' 


r 

J, 


nco 

+  rN  l20'(t)dt 


2a  r 


(t)  dt 


roo 

I 

*■  a2r^  j  p'(t)dt 


(3.  11) 


(3.  12) 

the  result 


0'(t)  dt 


(3.  13) 
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3 


The  above  result  may  ue  simplified  with  the  aid  of  the  following  identities. 


— —  f  t2<ft'{t)dt  =  i  + 


J 


2  +  x<p'(x)  -  <p  (x) 


(3.  Ha) 


1  2 

’°°  2  x 

-  I  te  <p'  (t)  dt  =  <p'  (x)  =  — — - - 

*/2it  J  >JZir 


1 


(3.  14b) 


-Lf  . 

\/2it  | 


1.2 


dt  =  <p  (x) 


^  0 


(3.  14c) 


The  resultant  simplification  is  the  following. 


2  ^  +  rN^  1 

<  6  >  - 2 -  ~  2  '  a  0' (a)  +  0  (Q)  +  2a^'  (a)  1  a* 

P 


■i  +  0(a) 


+  rN  |  1  *  a (°)  "  0  (°)  ~  2a  <p'  (a) 


+  [a2  i-  0(a)]  | 


(3  15) 


The  above  result  may  be  simplified  to  obtain  the  mean  square  outage  in  the 
following  form. 


(i 


<  0  >  =  p" 


rN  r  ' 


2(1  +  rN) 


r  “1 

,,  ,  2,  rN  2 

(1  la)  -  - - r-p  p 

N 


ja0'(i)  +  (1  +  a2)  <p( a)j  (3.  16) 


i  -6 


» 


where  a,  <p( a),  and  p'(a)  are  defined  in  Equations  (3.  12)  and  3.  14). 

/\  7*  2 

Mow  a  =  -  may  be  chosen  to  minimize  <0  >;  that  is,  choose  a  to 

prN 

satisfy  the  equation 


0  <e2> 

9  a 


0 


(3.  17) 


then 


(rN  +  «  =  <4  - 1} 


a  <j>"  (a)  +  (2  +  a^)  <p'  (a)  +  2 


a  <t>  (a) 


This  equation  may  be  simplified  by  noting  that 


p"  (a)  -  -  a  p'  (o.) 


(3.  19) 


Then  the  condition  that  the  mean  square  outage  is  a  minimum  is  expressed  by 
the  following  equation: 


■A'  (al  r  a  p(n) 


(3.  20) 


Tables  giving  the  normal  error  function  <p(a)  as  well  as  the  derivative  p'(a)  are 
given  in  the  Mathematical  Tables  from  the  Handbook  of  Chemistry  and  Physics. 


\ 


Tne  above  method  and  statistical  calculations  were  furnished  by  R.  A.  Moore. 
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When 


The  solution  to  Equation  (3.  20)  gives  a  as  a  function  of  r^. 

A  i* 

ft  “  -  is  known,  the  quantity  Ar  may  be  evaluated.  Then  the  loading  mixture 

prN 

ratio,  r,  ,  may  be  evaluated  as 


rL 


N 


Ar 


=  rN(l  -  up) 


(3.21) 


Minimum  mean  square  outage  loading  results  when  the  missile  stage  is  loaded 
on  the  basis  of  r^,  in  accordance  with  Equations  (3.  2),  (3.  20),  and  (3.  21). 

Table  V  gives  a  tabulation  of  Equation  (3.  16)  for  the  case  when 
r^,  =  2.  25.  Note  that  minimum  mean  square  outage  occurs  when  n  ~  0.  6.  The 
minimum  value  could  also  be  obtained  from  Equation  (3.  20). 

A  comparison  of  standard  loading  and  minimum  mean  square  outage 

loading,  when  r^  =  2.  25,  can  be  made  using  Table  V  Standard  loading 

requires  that  the  missile  be  loaded  with  n  =  0.  With  minimum  mean  square 

outage  loading,  a  -  0.6  is  used.  The  comparison  will  be  made  assuming 

<r ,  =  <r  =  0.5 (r„  -  0.25%.  Then  p  -  0.75%.  When  standard  loading 
1j  r  r 

is  used, 


\/<e^>  =  0.  536p 


=  0.  402% 


(3.  22) 


When  minimum  <9  >  loading  is  used. 


6J> 


=  0. 44  lp 


-  0.  331  Vo 


(3.  23) 
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SECTION  4 


PROBABILITY  THAT  THE  OUTAGE  IS  LESS  THAN  A  GIVEN  VALUE 

When  the  guaranteed  range  of  a  missile  stage  is  calculated,  it  is  assumed 
that  an  excess  of  propellant  is  on  board  to  allow  for  the  existence  of  outage. 

The  amount  of  additional  propellant  depends  on  the  probability  that  the  actual 
outage,  6,  will  be  less  than  some  reference  outage,  denoted  by  0q.  Recall  that 
both  0  and  0q  are  measured  as  percentages  of  the  total  propellant  loaded  per 
stage. 

It  is  shown  in  Equation  (2.  7)  that  the  fuel  and  oxidix_i  outage,  expressed 
as  a  fraction  of  total  stage  propellant,  is  given  by 


60„  ,  /6M„  ,  6  M. 

F  1  F  6  r  L  Ar 

M  '+r  M  +r~M  +"r 
•‘N  ‘  rN  \MFN  N  LN  N 


60L  -rN  /<>MF  .  6  r  *M 


L  ,  Ar 


MN  1  '  rM  \MFN  rM  MLN 


1  —  1 
'«) 


66v 

That  is,  and  ■..,■■■  -  ar.  i>«.  rtecFy  eoi-' since 

MN  mn 


60l  60f 
7^7  ‘  ’rN 


4-1 


Tlx-  actual  outage.  9,  may  then  be  written  in  terms  of  the  random  variable,  X, 
as 


0  = 


X 


1  +  r 


if  X  >  ft 


N 


"  ^  rM 

9  -  .. - -  if  X  <  9 


)  +  r 


N 


(4.4) 


where 


X 


6Mt 


M 


FN 


6Mr 


M 


LN 


(4.5) 


Tito  probability  that  0  <  9q  is  then  the  probability  that  X  occurs  along  the 
strip 


-  (l  +  rN) 
rN 


%  <  X  <  ('  '  «o 


(■I.  <•) 


Reference  to  the  figuie  below  will  help  rlai  ifv  Ihe  above  result. 
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Then 


P(6  <  0Q) 


+  rN>  60 


-  <*  +  V  eG 

rN  ' 


f  (x)  dx 


(4.  7) 


where  f  (x)  is  the  distribution  of  x.  Note  that  f  (x)  is  normal  with  nonzero 
mean.  That  is, 


f(x) 


_ 1__ 


2P 


(x 


Ar .  2 

rN 


(4.  8) 


where 


2 

P 


2 

°F 


(4.9) 


The  above  statistical  determination  was  furnished  by  R.  A.  Moore.  When 
Equation  (4.8)  is  substituted  into  Equation  (  '  71  with  a  change  L*  the  integration 
variable,  the  result  is 


P(9  <e0)  =-p 


f  < 1  +  rN) °o 


■l  +  0O 


N 


Ar  \2 
prN/ 


du 


(4.  10) 
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Equation  (4.  10)  may  be  written  in  the  form 


f 


p(e  <  e 


o' 


l 

n/2tt 


1  u 

2 


e  du  =  0(«q)  +  0(uj) 


where 


<p  (u)  = 


r 


1  2 
'  2  x 

e  dx 


U1  = 


{1  +  rN}  90  Ar 


P  r 


N 


(1  +  rN)_%  +  Ar 


0  i>  r 


N 


pr 


N 


For  example,  if  the  mixture  ratio  -  1,  Ar  0,  -  j  P  • 


ilien 


n 


<  l  o)  rr  -L 


p  (9  <  4  p) 


- 


1  2 

du  =  0. 6826  . 


(4.  11a) 


(4.  lib) 


(4.  11c) 


{4.  12) 


4-4 


I  hat  is,  in  this  cast  *he  actual  outage  will  be  equal  to  or  less  than  the  rms 
"ltage  about  two-thirds  o!'  the  time.  This  example  is  included  only  to  show  the 
inner  consistency  of  the  theory  in  an  obvious  case. 


Equation  (4.  11)  contains  the  parameter  a  = 


Ar 


pr 


used  in  Section  3.  The 


N 


parameter  a  is  still  adjustable  and  may  be  chosen  to  maximize  P  (9  <  9^). 

This  procedure  will  give  a  value  for  the  parameter  a  that  is  different  from  the 
value  obtained  when  the  mean  square  outage  was  minimized.  This  procedure 
will  then  lead  to  a  different  loading  mixture  ratio.  When  the  loading  is  based 
on  maximizing  P  (9  <  6^),  the  mean  square  oc'age  will  be  greater  than  with 
minimum  mean  square  outage  leading. 

It  is  clear  that  the  value  of  that  maximizes  P  (0  <  9_)  is  determined 

PrN  0 

by  equating  Equations  (4.  lib)  and  (4.  11c).  Then 


(1  + 


lN'"0 


_Ar 
o  r 


(1 


rN*0O 


Ar 


N 


P  r 


N 


P  r. 


N 


or 


\_r 

rN 


(4.  13) 


The  value  of  a  is  obtained  as 


u 


rN  ~  1 

2PrN 


(4.  14) 
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Substituting  the  calculated  value  of  a  into  Equations  (4.  lib)  and  (4.  11c)  gives 
the  result 


(1  +  9q 


VI,.  =  U  ,  =  - - 

0  1  ipr 


N 


Then,  the  maximum  value  of  P(8  <  0^)  is 


po 


FI6<V^Jo  ° 


1  2 
2  u 


du  =  Z<p  (u( 


where  u^  is  defined  in  Equation  (4.  15). 

The  mean  square  outage  that  results  when  P  (9  <  9^)  is 
be  calculated  using  Equations  (3.  16)  and  (4.  14). 

The  probability  that  the  outage  is  less  than  one  percent, 
be  compared  for  each  of  the  three  methods  of  loading.  As  a  i  c 
the  case  where  r^  -  2.25,  (r^,  -  0.25%,  o-j  0.5%  and  or 
before, 

f~T  ,  ?.  ,  2  . 

-N_  •„  +  <r-  +  r  =0.  75% 


Standard  loading: 

By  definition,  Ar  =  0,  and  from  Equation  (4.  11),  Uj  ~ 


(4.  la) 


,)  (4.  16) 


;  maximized  may 

P  (0  <  1%),  may 
va'nplo,  consider 
•  0.  5%.  Then,  as 


(4.  17) 


4.  33,  Uq  =  1 .92t>; 
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henco. 


P{0  <  1%)  = 


•m 


f4.  33 


>-  1.926 


1  2 
2  u 


du  =  0.  973 


(4.  18) 


Minimum  mean  square  outage  loading: 

When  =  2.25,  a  =  0.60;  hence,  from  Equations  (4.  11),  Uj  =  3.74 
and  Uq  =  2. 527,  and 


P(0  <  1%)  =  0.9941 


(4.  19) 


Maximum  P(9  <  1%)  loading: 

Substituting  r^  =  2.25,  9Q  =  1%,  p  -  0.  75%  into  Equation  (4.  15) 
gives  Uq  =  3.  13.  Then,  from  Equation  (4.  16), 

1  2 
-  1  u 

e  6  du  =  0.9983  .  (4.20) 


P(9  <  1%)  - 


2 

^~Jo 


The  mean  square  ou*aee  that  results  w  u-n  R{ «-  *  i"M  ?s  minimized  is 
easily  calculated  using  Equation  (3.  16)  (wmui  is  tabulated  in  Tab!  a  V  when 
r^  =  2.25).  From  Equation  (4.  14).  a  =  1.2;  hence,  the  resuitanr  rms  outage 


'  Similarly  it  may  be  shown  that 

Pi  9  ■-  C.  49?06)  ■■  9.  80 
P(9  <  0.  526%)  =  0.  90 
P(0  <  0  626%)  =  0.  95 


P(9  <  0.  744%)  -  0.  98 
P(9  <  0.822%)  •:  0.99 
P(9  <  1.  046%  =  0.  999 
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.  /  2 

i<;\/<  6  >  -  0.  375%.  It  will  be  noted  from  Equations  (3.  22)  and  {3.  23)  that  the 
ti'S  outage  is  increased  over  the  value  obtained  foi  minimum  mean  square 
outage  loading  but  less  than  that  obtained  with  standard  loading. 

The  loading  mixture  ratio  is  again  computed  using  the  result: 


M 


0  -  ^) 

N 


(4.21) 


When  P(9  <  9^)  is  maximized. 


^  r 

-  is  given  by  Equation  (4  13),  that  is, 

rN 


r 


L 


r 


N 


1) 


(4.  22) 
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SECTION  5 


GENERAL  DISTRIBUTION  AND  STATISTICS  OF  OUTAGE 


In  general,  the  outage  frequency  distribution,  f(4),  may  be  obtained  by 
differentiating  Equation  (4.  11)  as 


Hence, 


m  = 


3P(9  <  e0) 


8  9, 


J  %  *  9 


(5.1) 


m  = 


\l  Ztt 


I  2  I  2 

'  2  U1  (1  +  r  )  ' 2  U0  (1  +  r  ) 

e  -  i  e  11 


pr 
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1  +  r 
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'JZtt 


<  e 


2p 


(1  +  rN)  e  -  — 
N  rN 


'N 


-22 

2p  rN 


(1  +  rN)  0  +  Ar 


.1-  c 
1* 


{5.2) 


Mean  outage: 

When  the  distribution  of  outage  is  known,  the  mean  outage,  9^,  is 
obtained  as 


poo 

&N  =  j  0f(9)d0 


(5.3) 
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When  f(0)  from  Equation  (5.2)  is  substituted  into  Equation  (5.3), 

2 


1  x  r 


e 


N 


N 


p  'n/  2rr 


17 


,22 
1  „  2p  rN 


I- 


(1  +  rN)  0  +  Ar 
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+  — t-  e 
rN 


The  result  of  the  integration  is 
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P 


{• 


'(a)  +  a  $i(u) 


where 


-  At- 
"  prN 
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du 


*'(  a)  =  ■— 

\>  2ir 


dO 

J  (5.4) 


(5.5) 


(5.6) 


5- 


A  loading  mixture  ratio,  ,  may  !><■  obtained  in  such  a  way  that  9^  is  a 
minimum;  then  Ar  will  jc  a  solution  of  the  equation 


MN 

oTatt 


The  result  of  differentiating  Equation  (5.5)  is  obtained  as 


where 


0"(a)  4-  a ^'{a)  +  0(o) 


2T7~rn 


1 

-  a  e 

2it 


1  2 
2U 


(5.7) 


(5.8) 


(5.9) 


Since 


^"(<> )  -  -  a  V>'(<i ) 


(5.  10) 


it  follows  that 


0(ft) 


_1 

~> 


r 


N 


1 


(5.  11) 


t) 


specifies  the  value  of  a  =  — ■ —  that  will  minimize  the  mean  outage,  0N< 

^rN 

When  the  loading  is  chosen  to  minimize  the  mean  outage,  Equation  (5.  11) 
will  apply.  In  this  case,  a  simplification  in  Equation  (5.5)  is  obtained.  The 
result  is 

»(VJ  -  P0'(«)  (5.  12) 
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:c  variance: 


The  outage  variance  (about  the  mean),  cr^,  is  calculated  as  follows: 


The  definition 


(5.13) 


eN  =  <  0  >  (5.  14) 

was  used  in  the  above  derivation. 

2 

The  quantity  <0  >  is  the  mean  square  outage  derived  in  Section  3. 

This  same  expression  may  be  alternately  obtained  as 


f  oo 

I  e2  f(9)  de  (5. 15) 

Jo 


When  the  outage  frequency  distribution,  f(0),  from  Equation  (5.2)  is  substituted 
into  Equation  (5.  15)  and  the  integration  is  performed,  the  result  is 


(5.16) 


7 

<  0"  > 


r2  i 
N  ' 


(1  +  a2) 


2(1  +  r  N)' 


I  7 

j  a  t  (l  +  a  ')  ip(a) 

t_ 

Note  that  Equation  (5.  16)  is  in  agreement  with  Equation  (3.  16).  The  expression 
for  mean  square  outage  is  simplified  when  the  missile  is  loaded  30  as  to 
minimize  the  mean  outage.  In  this  case 


<  92  >  =  p2 


(1  +a  )r 
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(i  +  rN' 


r.T  -  1 

y  ■  -+-[■  »  *'(<*) 
rN 


(5.17) 


The  outage  variance  when  the  missile  is  loaded  to  minimize  the  mean  outage  is 


o' 
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*  qS>N 

(1  +  r  )2 

l'  n' 


rN  -  1 
rN  +  1 


a  ^'(u) 


a 


(5.18) 


When  the  missile  is  loaded  in  any  other  way  (i.  e. ,  such  that  Equation  (5.  12)  is 
2 

not  satisfied),  then  <Tg  may  be  obtained  using  Equation  (5.13).  That  is, 


2  . 


<  e2  > 


e 


2 

N 


z  gi "on  by  Equals  (5.  16)  and  0^  is  given  by  Equation  (5,  5). 
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where 


Pro',  '■■bility  that  the  outage  is  fuel  or  oxidizer 

The  probability  that  the  outage  is  either  fuel  or  oxidizer  may  be 
h-termined  from  Equation  (4.4).  When  X  >  0,  the  outage  is  fuei;  when  X  <  0, 
fhe  outage  is  oxidizer.  Mow  X  is  a  random  variable  with  mean, 


<  X  > 


Ar 


The  distribution  of  X  is  given  in  Equation  (4.8).  The  probability  that  the  outage 
is  fuel  is  equal  to  the  probability  that  X  >  0.  Then: 


P(X  >  0)  = 


ra>  .  i_  x  -  *L) 

1  /  -  2p2  V  rN/  dx 


«/2iF 


*•  •'0 
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^  r 
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X 

—  r 


au 


(5.  19) 


where  <p(a)  is  defined  in  Equation  (5.  f>).  The  outage  can  only  be  fuel  or 
oxidizer.  That  is, 


P(X  t;  ■  )  4  P{X  <  0)  =  1 
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or 


P(X  <  0)  =  I  -  t{a)  .  (5.20) 

Hence, 

P{0  =  fuel)  =  —  +  0(a) 

P(0  =  oxidizer)  =  —  -  0(a)  .  (5.21) 


Note  that  when  the  ratio  of  the  two  probabilities  is  ^qual  to  r^, 


j  +  <t>M 

I  '  *( a) 
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N 


or 


0(a)  = 


_1 

2 


rN  ~  1 
rH  +  1 


(5.22) 


Equation  (5.22)  is  identical  with  Equation  (5.  11),  hence  ihe  above  assumptions 
pertain  to  minimum  mean  outage  loading. 
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SECTION  6 


BIASING  TO  MAXIMIZE  MISSILE  PERFORMANCE 

A  particular  missile  configuration  is  usually  required  to  place  a  payloaci 
into  a  ballistic  orbit.  A  single  parameter,  V,  may  bv  examined  to  determine 
whether  or  not  the  given  missile  configuration  will  have  the  capability  of 
achieving  the  mission.  The  discussion  in  this  section  will  assume  that  V 
refers  to  the  burnout  velocity  magnitude.  No  loss  in  generality  is  introduced 
by  this  assumption  because  V  could  equally  wcl'  refer  to  the  burnout  energy, 
the  total  range  at  earth  impact,  or  any  other  quantity  with  the  desired  property. 

The  nominal  velocity  of  the  missile  at  burnout,  V,^,  is  a  function  of  many 
variables,  for  example 


VN  =  VN  *MN’  FN’  IN’ 


•  V 


(6.1) 


where 


M^  =  expected  liftoff  weight 
F^  =  expected  thrust 
Ir_  -  expected  specific  impulse 
9^  =  expected  mean  outage 

The  three  dots  in  expression  (6.  1)  refer  to  other  variations  not  specifically 
listed.  These  include  trajectory  shaping,  missile  steering  with  guidance, 


Expression  (6.  H  may  be  expanded  in  a  Taylor's  series  to  obtain 


V  =  VN  +  A  AM  +  B  AF  +  C  AI  +  •  • .  -  D  A0  (6.2) 


where 


AM  =  liftoff  weight  variation 
AF  =  thrust  variation 
AI  =  specific  impulse  variation 
A0  =  0  -  0^  =  outage  variation 


and  A,  B,  C,  •  •  •  ,  D  are  the  corresponding  partial  derivatives.  That  is, 


A  = 


9V 

9M 


(6.3) 


where  the  partial  derivatives  are  evaluated  at  burnout.  Note  that  the  sign  on  D 
is  chosen  in  expression  (6,3)  such  chat  D  >  0.  It  is  assumed  that  AM,  AF, 

AI,  >-•  ,  A0  have  zero  mean  but  not  necessarily  that  these  variations  are 
normally  distributed.  Furthermore,  it  is  assumed  that  M^,.  F^,  1^  +  . . .  are 

specified  by  the  missile  cc—«';~u..tri.<j«.  This  is  not  the  case  for  0^.  The  value 

for  0^  depends  on  the  missile  loading  procedure  (or  how  the  PU  Bystem  is 

biased). 

Equation  (6.2)  may  be  simplified  by  defining  a  nominal  final  velocity, 

* 

V^,  that  does  not  depend  on  0^,  as 


+  D9N 


(6.4) 
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The  quantity  V^  is  v/'tained  by  using  the  nominal  expected  values  for  M,  F, 

I,  •  •  •  ,  and  allowing  the  missile  to  burn  out  with  zero  outage.  Then 

V  =  V*  -  D  eN  t  A  AM  +  B  AF  +  C  AI  +  •  •  •  -  D(9  -  eN)  (6. 5) 

It  is  clear  that  if  the  performance  of  the  nominal  missile  is  to  be  maximized, 
the  missile  should  be  loaded  to  maximize  9^.  The  loading  mixture  ratio, 

r^  =  r^  -  Ar,  that  will  maximize  0^  is  obtained  in  Equation  (5.11). 

It  may  be  desired  to  maximize  the  probability  that  th"  missile  performance 
(as  measured  by  the  parameter  V)  is  greater  then  some  fixed  value,  Vq.  The 
development  in  the  following  paragraphs  applies  to  this  case. 

It  is  desired  that  the  missile  achieve  a  required  velocity,  Vq,  with  a  high 
probability.  The  probability  that  V  will  be  greater  then  Vq,  F(V  >  Vq),  may 
then  be  written  as 


P(V  >  Vq)  = 


(6.6) 


where  f  (V)  is  the  frequency  distribution  of  V.  It  is  assumed  there  are 
sufficient  random  variables  ,ua!  I>ir,  (6  5)  to  .justify  the  use  of  the  central 
limit  theorem.  Then  the  distribution  of  V  will  be  normal  .*  1th  m"*" 


<V> 


D  9 


N 


(6.7) 
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and  variance 


<r 


2 

V 


( 


Ta  AM  +  B  AF  +  C  AI  + 
L 

2  ,  _  2  2  ,  nZ  Z  , 

+  B  ^f  +  C  ^  I  + 

<r  2  +  B2tr2  +  C2crZ  + 

M  5  I 


-  D  (9 


+  D2<r  2 


(6.8) 


where  cr^,  cr^,,  cr^,  •  •  •  ,  cr^  denote  the  one-sigma  variations  about  the 

expected  mean  of  the  variables  M,  F,  I,  •  • •  ,  9.  Note  that  all  cross 
correlations  were  assumed  to  vanish.  Hence,  Equation  (6.6)  may  be  written 
as 


P(V  >  VQ)  = 


(6.9) 


When  the  substitution 


u 


(6. 10) 


is  made  in  Equation  (6.  9),  the  following  result  is  obtained 


'<  / 


V  = 


(6.11) 
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where 


u 


0 


VN  +  D9N 


(6.12) 


It  is  desired  that  the  missile  stage  be  loaded  so  that  P(V  >  Vq)  is  a  maximum. 

It  is  clear  from  Equation  (6.  11)  that  this  is  equivalent  to  biasing  such  that  Uq 

9uq 

is  a  minimum.  When  — —  is  formed  and  set  equal  to  zero,  the  result  is: 

oAr 


!!n 

9Ar 


V*  +  D9,.  ocr 

_ IN  _ V 

D<Ty  9Ar 


(6. 13) 


When  Ar  is  obtained  as  a  solution  of  Equation  (6.  13),  the  loading  mixture 
ratio  r^  =  r^  -  Ar  is  obtained  that  will  maximize  P(V  >  Vq). 

Equation  (6.  13)  may  be  written  in  several  alternate  forms.  When  it  is 
assumed  that  AM,  AF,  A  I,  etc.,  are  independent  of  Ar,  the  following  result 
is  obtained: 


3<r 

V 

°V  OAr 


i  0rl 


2  9Ar 


Hence,  Equation  (6.  13)  may  be  written  as 


99.. 

N 

9Ar 


'N 


D°N  S 


(0) 


9Ar 


(6.  14) 


(6.15) 
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Now 


<  e2  > 


(6. 16) 


Hence,  Equation  {6.  15)  becomes 


09.. 

N 

9Ar 


(1  + 


2*  V  = 


<  e2  > 

9Ar 


where 


X  = 


+  D6N 
O) 


(6.1?) 


(6. 18) 


Equation  (6. .17)  is  best  solved  for  Ar  by  iterative  methods.  The  resulting 
value  of  Ar  will  specify  the  loading  procedure  that  maximizes  P(V>VQ).  When 

2  2 

Ar  is  known,  the  quantities  G^,  <0  >,  and  cr^  may  be  calculated.  Then 
2 

cr y  (9)  and  P(V  >  Vq)  may  be  evaluated. 
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SECTION  7 


l 


GENERAL  DISCUSSION  OF  A  PROPELLANT  UTILIZATION  SYSTEM 


> 


Outage  is  minimized  in  the  Thor  and  Titan  missiles  using  an  open-loop 
system.  This  method  requires  that  close  tolerances  be  maintained  in  propellant 
loading  and  the  burning  mixture  ratio.  The  purpose  of  a  PU  system  (as  used 
in  the  Atlas)  is  to  cause  the  fuel  and  oxidizer  to  be  expended  simultaneously 
using  closed-loop  measurement  and  control. 

The  operation  of  a  PU  system  requires  th  it  the  remaining  amounts  of 
fuel  and  oxidizer,  F  and  L,  respectively,  be  measured.  Then 

F*  =  F  -  f  (7.  1) 

L*  =  L  -  i  (7.2) 


where 


F,  L  = 
f.  t  = 


remaining  fuel  and  oxidizer  measured 

actual  fuel  and  oxidizer  remaining 

errors  in  the  fuel  and  oxidizer  measurements 


The  measurement  errors  f  and  l  are  assumed  in  this  analysis  to  be 
constants  for  any  one  flight.  The  units  are  mass  (pounds)  of  fuel  or  oxidizer, 
respectively.  The  frequency  distribution  of  these  errors  for  an  ensemble  of 
flights  is  assumed  to  be  gaussian  with  mean  zero.  The  measured  tank  mixture 
ratio  may  then  be  defined  as 


V  Jj 

R  =-* 
F 


(7.3) 
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If  the  actual  burning  mixture  ratio,  R,  is  controlled  to  equal  R  ,  then  the 
equation  of  operation  of  the  PU  system  may  be  written  as 


$ 

R  =  R  +  6r 


(7.4, 


where 


6r  =  gaussian  constant  error  in  controlling  the  mixture  ratio 
L  =  actual  oxidizer  flow  rate 
F  =  actual  fuel  flow  rate 

Equation  (7.  4)  may  be  written  as 


F  F  -  f 


It  may  now  be  assumed,  without  loss  in  generality,  that 

T  "  F. a  constant 

t\i 


(7.5) 


(7.6) 


or 


F 


F 

0 


fn1 


(7.7) 


where 


actual  initial  mass  of  fuel  loaded. 
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Equation  (7.5)  may  t«.  a  be  integrated  to  obtain  the  functions  L(t)  and  R(t). 
The  result  of  this  integration  shows  that  R(t)  varies  as  a  function  of  loading, 
measurement,  and  control  errors. 

The  above  unpredictable  variation  in  R(t)  can  be  avoided  by  controlling 
the  burning  mixture  ratio,  R(t),  to  the  nominal  burning  mixture  ratio,  r^. 

In  this  system,  the  flow  rate  of  the  oxidizer  or  the  fuel  is  controlled  to  be  pro¬ 
portional  to  the  difference  between  the  measured  tank  mixture  ratio  and  the 
nominal  burning  mixture.  The  operation  of  this  system  is  described  by  the 
equation 


L  -  LN  =  K  (R*  +  6r  -  rN) 


(7.8) 


where 


K  =  constant  of  proportionality 
■Lj^,  =  nominal  oxidizer  flow  rate 

LN 

rM  =  nominal  burning  mixture  ratio  -  - — 

F 

XN 


Note  that  Equation  (7.8)  is  written  assuming  that  the  oxidizer  flow  rate  is 
adjusted.  The  fuel  flow  rate  may  be  chosen  to  be  nominal  without  loss  in 
generality.  Note  also  that  if  K  =  -  F^,  the  system  described  by  Equation  (7.  5) 
result?.  Since  this  system  may  be  shown  to  have  some  undesirable  character¬ 
istic.  ,  tu°  constant  K  will  be  chosen  such  that 


K  /■-  F 


N 


(7.9) 
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In  general,  K  may  be  written  as 


K  =  (xFn,  (i  P  I  (7.10) 

where  p  may  be  referred  to  as  the  PU  system  gain. 

The  system  operation  described  by  Equation  (7.8)  closely  approximates 
the  Convair  PU  system  used  on  the  Atlas  missile.  Equation  (7.8)  may  be 
written  as 


1  -  LN  =  (?-TT+  6r  -  Tn) 

or 

L  -  -JL  L  =  LN  -  (ii'N  (rN  -  6r  +  — —  j 

Equation  (7.  12)  may  be  integrated  to  obtain  L(t)  when 


(7.11) 


(7.12) 


F 


F0+  FNt 


(7.H) 


and  Fj^,  is  a  constant.  The  result  of  this  integration,  when  n  M,  is  given 
below. 


(7.14) 
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where 


£ 


f> 


Lp  =  actual  initial  mass  of  oxidizer  loaded.  (7.  15) 

Equations  (7.  13)  and  (7.  K)  describe  the  PU  system  operation  under  the 
assumption  that  the  fuel  flow  rate  is  a  constant.  Now  consider  that  the  PU 
system  has  been  operating  for  some  time.  Then 

F  <  <  Fq  (7.16) 


and 


L  <  < 


(7.17) 


4 


Hence,  if  the  PU  system  gain,  p,  is  large,  then  e  will  be  small  compared 
with  the  first  term  in  Equation  (7.  14).  The  magnitude  of  e  is  further  reduced 
if  an  attempt  is  made  to  load  the  propellant  in  accordance  with  the  nominal 
burning  mixture  ratio  r^.  Then,  since  c  is  a  second-order  term, 

L  - '  ■  ('h  -jrrr 6r) |F  -  *'  •  i7-18) 

Satisfactory  operation  of  the  PU  system  will  be  obtained  for  any  large  value 
of  p.  For  example,  p  may  be  chosen  in  the  range 

2  <  p  <  5 
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Equation  (7.  18)  was  derived  from  Equation  (7.  11),  which  requires  the  oxidizer 
i!ow  rate  to  change  linearly  with  the  error  signal 


In  any  practical  system  there  is  a  limit  to  the  allo%vable  change  in  the 
oxidizer  flow  rate.  Hence,  in  the  initial-operation  phase  of  the  PU  system,  the 
oxidizer  control  valve  may  be  hard  over  (one  way  or  the  other)  and  Equation 
(7.  11)  will  not  apply.  It  is  assumed,  however,  that  the  tank  mixture  ratio  will 
eventually  be  controlled  so  that  Equation  (7.  11)  does  apply.  If  the  PU  system 
gain  is  not  too  high,  the  system  should  then  continue  to  h"  linear.  If  p  is  too 
high,  an  undesirable  limit  cycle  may  result  because  of  inertia  in  the  valve 
hardware. 

The  quantities  F^  and  do  not  appear  in  Equation  (7.  18).  It  then  may 
bo  concluded  that  outage  is  not  a  function  of  initial  loading  errors.  As  the  pro¬ 
pellant  approaches  depletion,  the  mixture  ratio  control  error,  6r,  becomes  a 
second-order  error,  and  (independent  of  the  PU  system  gain)  Equation  (7.  18) 
reduces  to 


L  -  i  =  rN  (F  -  f)  .  (7.20) 

Note  that  when  F  and  L  have  been  sufficiently  reduced.  Equation  (7.20)  re¬ 
quires  that 


Because  of  the  fast  reaction  lime  of  the  system  (with  high  p),  this  equation  will 

be  IrU’’  eve:,  wi.en  ,  and  F  h..\e  small  variations  with  time. 

1\ 
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The  propellant  measurement  system  may  be  designed  such  that  L  and  F 
become  ambiguous  when  only  a  small  amount  of  propellant  remains  in  the  tanks. 
If  this  is  the  case,  it  will  he  assumed  that  the  system  goes  open  loop  prior  to 
tne  occurrence  of  the  ambiguous  measurements  with  R  =  r^  until  one  of  the 
propellants  is  exhausted.  A  small  control  error  during  this  process  may  be 
neglected  because  only  a  small  amount  of  propellant  remains.  When  this  pro¬ 
cedure  is  followed,  no  change  in  Equation  (7.20)  is  required. 

The  outage,  whether  it  is  oxidizer  or  fuei,  may  now  be  determined  from 
Equation  (7.20)  as  a  function  of  the  measurement  errors  f.  and  f  which  were 
assumed  to  be  constant  in  this  equation.  Because  of  the  fast  reaction  time  of 
the  system,  f  and  i  may  be  time-varying  quantities.  The  values  of  f  and  t  to 
be  used  in  calculating  outage  are  the  values  of  f  and  f  that  persist  when  the 
propellant  is  nearly  expended. 
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SECTION  8 


OPERATION  OF  A  PROPELLANT  UTILIZATION  SYSTEM 

Equation  (7.  20)  may  be  used  to  compute  the  outage  that  results  when  no 
biasing  technique  is  employed.  When  the  fuel,  F,  is  exhausted,  the  outage  is 
oxidizer  and 


0L  =  i  -  rN  f  >  0  (8.  la) 

When  L  =  0,  Equation  (7.20)  may  be  solved  for  the  fuel  outage,  0  ,  with  the 
result 


0 


F  = 


(P.lb) 


Equations  (8.  la)  and  (8.  lb)  may  now  be  compared  with  Equations  (2.  9a)  and 

(2.9b).  The  mean  square  outage  may  again  be  calculated  using  Equations  (8.  la) 

or  (8.  lb)  whichever  gives  a  positive  result.  It  will  be  assumed  that  f  and  i  are 

uncorrelated  and  normally  distributed  with  mean  zero.  Then  the  outage  is 

2, 

equally  likely  to  be  ail  fuel  or  all  oxidizer.  The  mean  square  outage,  <0  >, 

written  as  a  fraction  of  the  total  propellant,  may  then  be  evaluated  as 


2 

<  9  > 


Z2  4-  z2 
Zf  .  zi 


(8.2) 


2  2 

where  Z.  and  Z^  are  mean  square  fuel  and  oxidizer  outages  expressed  as 
fru.tio'  s  et  the  t<.Uu  propel'ant.  These  quantities  may  be  calculated  from 
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$ 


q  a.ions  (8.  1)  by  so  >  .ring,  averaging,  and  dividing  by  the  square  of  the 
nominal  mass  of  propellant,  as: 


Z 


2 

f 


and 


where 


(1 


+  rN> 


K  +  # 


<4) 

(1  *  rN)2 


//_±_ 

\j  mln 


(8.3) 


(8.4) 


(8.5) 


(8.6) 
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and  M,  ,  ;ind  ivl .  ..  again  denote  the  nominal  total  mass  of  propellant, 

N  l*  N  LN  n 

fuel,  and  oxidizer  ioaded.  The  symbol  <  •■•  >  indicates  an  ensemble  average 
over  a  large  number  of  flights.  When  the  above  results  are  substituted  into 
Equation  (8.2),  the  expression  for  mean  square  outage  is  obtained  as 


T 

<  0“  > 


rN  +  1 


2d  +  rN)‘ 


2 

q 


(8.7) 


where 


2  2,2 

q  <rf  +  crjf 

o*£  =  rms  fuel  measurement  error  (l<r),  expressed 
as  a  fraction  of 

(r^  =  rms  oxidizer  measurement  error  ( 1  cr) ,  expressed 

as  a  fraction  of  Mj^  (8.  8) 

The  expression  for  mean  square  outage  given  by  Equation  (8.  7)  assumes  no 
biasing  of  the  PU  system.  V. i.:u  is  the  ..y:  rom  operation  is  specified  by 
Equation  (7.  1 1). 

The  similarity  between  Equations  {8.  7)  and  (2.  14)  should  be  noted.  The 
equations  are  identical  except  that  p  has  been  replaced  by  q. 


The  operation  of  a  PIJ  system  can  be  improved  by  introducing  a  bias  into 
the  rnea> uronv. r:v  *.  cither  thi  fuel  or  the  oxidizer.  When  such  a  bias,  Af,  is 


introduced  in  the  lucl  measurement,  the  oxidizer  or  fuel  outage  expressed  as 
fractions  of  the  tot.  1  propellant  becomes 


rN 


1  4-  r 


N 


l  +  Af 

MLN  mfn, 


(8.9a) 


or 


1 


1  +  r 


N 


M 


FN 


l  Af 

mln  MFNj 


(8.9b) 


These  equations  should  be  compared  with  Equations  (4.  1)  and  (4.?.).  An 
analysis  nearly  identical  with  that  used  in  Section  4  will  derive  the  expression 
for  P(0  <  0q)  when  a  FU  system  is  employed.  The  result  of  this  derivation  is 


P(0  <  80)  = 


1 

^Zk 


(8.  10a) 


where 


V1  * 


(]  +  rN)  e0  Af 


qMT 


(8.  10b) 


V0  * 


’  rN^  90  +  Af 


qr 


M 


q  M. 


FN 


(8.  10c) 
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Equations  (8.  iO)  should  be  compared  with  Equation  (4.  II)  and  the  similarity 
noted.  These  equate  .is  will  be  the  same  if  the  following  identification  is  made: 

p  =  q  (8.1  la) 


Ar  _  Af 
PrN  c1Mpn 


(8.  lib) 


The  above  identification  may  bo  used  in  Equations  (3.  16),  (4.  11),  (4.  16),  (5.5), 
(5.  17),  etc.  Then  ail  of  the  relations  (derived  on  the  assumption  that  conven¬ 
tional  loading  techniques  would  be  employed)  may  be  extended  to  apply  to  PU 
system  operation.  Further,  it  should  be  observed  that  the  PU  system  fuel 
measurement  bias,  Af,  is  equivalent  to  the  fuel  loading  bias,  AMp,  discussed 
in  Equation  (3.  7a). 


Each  cf  the  loading  methods  described  in  Sections  3,  4,  and  5  result  in  a 

Ar 

value  for  the  parameter  a  = 


pr 


,  When  a  and  p(=  q)  are  known,  Ar  or 


N 


AMp  (=  Af)  may  be  determined  using  Equation  (8.  1  lb). 


Proper  operation  of  the  propellant  utilization  system  requires  that  the 

$ 

quantity  Af  be  subtracted  from  the  fuel  measurement,  F  .  The  oxidizer  flow 
rate  is  then  coni  rolled  according  lo  the  equation 


L  ‘  LN  “  1111 N 


Af 


+  6r  -  rN 


(8.12) 


rather  than  Equation  (7.  11). 

A.,  oxi  d .uiitrcmnii*  bia.?  Af,  could  have  been  determined  rather 
lhan  Af.  This  quantity  would  he  obtained  from  Af  as 
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Af 


(8.13) 


A  l  _ 

MLN  mfn 


In  this  case  the  oxidizer  flow  rate  should  be  controlled  by  the  equation 


LN  =  ‘tF 


N 


A  i 


+  6r 


’  N 


(8.  14) 


rather  than  Equation  (7.  11). 


8-6 


1 

2  Po  ^  £ 
4  o  H  « 

*g_-2  I 
S§«=  c 

a:--  J 
§  J  >-T 
Z.  <  <3  c» 

So-gs 

3gi^ 

c  .  *r 
„  £  o  «o 

^  <  o  v  rr 

r-  ~  O 
Zm>2  l*fL 
®< 


aS  U  *  •£  P 
1HQ  2i 


c  c  ~  „ 

:  C  v  ft  JL  «.c^  . 

U  S'.  5  E  e  o  ^  ~  .c  ® 

vi ‘•s®  “V 

i  S.of'Bg/^ef 
:s~  S~°  ?£,  *v  2  o 

*  v  «n  J5  J;  *3  .flu 

,  i">  w  C  Q,  0<  .  1  y  4>  O 

;  2-£p  £  “  =  «  3  2  a.? 
;Sg  rS-S23*’^ 
n»-jS  3|5g5- 
|||M»53ts|a| 

>•££.£  SflS  .c?S 

:s°11?s || S  s  § 

:  I  2^  Sa2^»“ 

JC'-CflC-"  -Sb,. 

5  «  .S  0  ©  —  -C  ••  O  ^ 

:  g  S  g  5  5  «*  &  c  y£  - 

S  S  *-N«  ~  <D  M 

V  o$:k<s*j^  ~  - 

<  n  “  O  u  .  <"5  O  V  ft  *3  O 
.-S.2  V  .ft  w  >  M>  C/ 

:  A  -  o  o  .  t>  e  u  **  w  ** 

'■-rt^ioc^3513 
*^a  3i^S^Tj  o  o  o 


^  k,  «i 

--  c 

&S  -  ~  3 

7  » 

^  J  S^f 

:35S-~ 

^2  O'? 

4« 

U  P  ^5  1 

^  ^  £  O  3 

*<8** 
M  Z  „  © 

2  M  x  S  tu 
^ o.<I 

j2.-*  t 

j  £  *  -  « 
IStj’S 

3  O  •-  S  g 

*  Id  •-.  -O  ° 

i.t-Q  22 


i>  5  “J  i 

5  |2S  1 

X  =  4,  -O  I 
p  -O  ~  £  - 
c  v  *»  E 
t<  >*«  SO 
o  ®  S  S' 
a  - ! 

§  ej;  - 

■2  -£  r; 

2dT  .  »• 
2  £  •  «  " 
tjr; 

~  c  rj  M 
o-E^c* 

o  o  *-5  ‘ 
*  y  2  g  = 
?  «  a.2  ; 

■*  y  4 

Z  g.jS'S  ■ 

t  c  J  o 
^  x  .2  *2  < 
m  «  c  c  • 
o  «  u  y  • 
ft**  0  >  { 
t>  «  -  c  ; 

\s s  8  * 

ns  • 
x  ,2  *5  »  ? 

H  5  jQ  3  ** 


*  »  *  E  c 

S^c  o  J3  ^  j: 

v  *»  u  r\  o  * 

JfS-s|2£ 

S.2|<-|v 

Si  2-5  .-  s  , 

e  5  §  «  3  ; 

5  So  g.«t 

r53  M-'l 

ti  <  c:  «  .  4 
—  .-o^i 

■-  .  -3  ^5  ft  c 


o«u  S  -  “*£ 

ir  **  tfi  «i  s 

0  o  «  S  * 
a.  %  c.2  v  i 

S.-St^T 

'J  m  ^  c£  c 
J  S3  N  fl.Sc 

S  <  2  3  -3 

S  .  rt  o  v  t 
.  x.  w  >  6 

>»  «  W  Od  -  fl 

4>  »  O  fl  “  *< 


Oq  ^ 

IfiPS4 

«fjv>^. 
w  <  W  ci  c 
So  /s  j  T 
-io  ; 

2>S<0- 


A  —  W  W  /  H  **  H 
,  f  c.  *’ ,  gf  .1 

•;  d  v  -  H  **  f  f4 1 

■  ■  .  -  H 

;  «52  E4vi 

jsss-! 

i;  “  f  **  ?i  «  fl  «  7' 

f.  “  <J  **  'A  O  *  ft. 

r2"5,!l  n. 
s  •  s  %  *'  v  i;  o  ji  ; 
■*  .  «  V  ,,  ^  xj  i 

,300^, 

:'Si2;..‘  =  5 

5fl'* 

-~  =0  S’-1 


2 _z : 

:  S  *  3  a  I  f  a 

jlihff? 

v  ^  i* 

:  ■  j7  *'  'J  ►  H 
l  £  J\  ft  «  ^ 


:9  - 

5  2.--  o 

"MO  ~  K 

'  ,  1 

\  -J  >-*r 

i  <  «c 

4  *J  ifi 

(§>i? 
i°  .gt 

:h*  ,‘r 

iis1'.4 


o  "  0  .  1  .'  f 

‘•ji.ts'j'S; 

go5^"5>i3< 

f^*"— "O'*  fiOc1* 
£»S 

0^  H  I- C «S> 

..  a  v  5  c  5  y  »  , 

i>  SiS^I  s  2  Sj 

£2r  .“Self; 

COjit  .  j  d  i 

H.gSdSS'-'S 

®  o  «"3  S  2  «  c  • 

*  y  r  ^  SI  ^  ^  v  ’ 

sr^s.oti) 

sli|s.§Ss; 

’Jf  «.2  u  «  2_j 

s*6gs!«5« 
s  I  §*5sli 

2'm  v  >.  -it .'  •• 

s  n  s:  £  t»  .oft* 

(-dJ  5-5*223 


<4  3  — 

?  C  " 
ft  w 

-  v  rt  . 

Va  4i  tft 

0  ~  dP 


Je  \  ^  JV 

? 

V  o  o 

3  Sf  a^. 

sixi 


U  *J  v  *» 
U  3  U  3 
*3  O  fC  ft 


2  i 

•V  J 

i  <  « <> 

<  5  cJ  c* 

^  P  Q  **> 

£  <>3 
=  <8*~ 
i^yjs  I* fj. 

1o5 

J  —  O  ol  y 
a  g  >  ~  2 

ir-j  »  t 

W  S  -2  n 

iHQ  T:Z. 


5  =  t;  v 

»  o  —  £  0  “  ■ 

IV  H  C-C  O  ^  v?  /:  ® 
5-5Wwi,'30>_*xy 

«  «  a  G  C  A  V®  ^ 

"■oirc'g^-2'®  =  a> 
=  3  s:  .2  £<2  «v|  .£. 
:  s&  e-J  ■  3  -j-g  S 
cScuSsSsS.'; 
g  =5*5  »£’»•= 

ii3*Pli^ 

pi -1 1 

f=S5lng--|*” 

S-Slo2i|.£”| 

:  2  £  1  "  «  “  JIT 

»  O  a-1  5-0*'  •  £ 

»  -2  v  w  .2  — •  *c  ■*  O  £ 

£  f,  32  «  in  z=” 

2  C  §  ^  S  3  «3  .  u  "* 
j  o  J  .  4  o  »  v  c  a 
2  u  .  je  "■*  >  ce  >  u 

-  *j  2 

33**'— .£.*-.’^5Or50 


a 

v 

>.*0 


<  r'  a2  £ 

2.  Hvo  - 

J  CJ  « 
•*  .2 
*  U  i!  c 
2“2=  = 


<  _ cc  S' 

Jg  Df 

rH£0-3 


•4  Z  „  o 

iMi'x 

i2.-lin 

jy.*~  « 
Ut'J 
)0  -»  c 
*u  ^3  0 

liHQiiv 


M  «  ^  ,  5 

=  y  o  w  —  .  -  o-  - 

2  S  *»T  p  v  SL  *  *■  «  J*  a 

S  IS  S'£  2J  oi-  i 
c 

-  4,-0  «  c  V  0 

a  ~  5^r-.2  Fu.J2v«  i 

!»  «  3  w  U  *j  ft  ; 

ca  a  3  .  u  .' ©  - 

2  e  2~ ’C  c  *  n  • 

S-S'^«3S5«a: 

sSsSoS^I 

£  *5  <  C  £  t  *  •  * 

c  j  .;«  •-  o  #  >-  u  i 

Oe^^-O.SJ- 

£  rt  S*  $  *  S'"?  ’ 
c  JvS  c  <e  .  '-  .  w  *  1 
0  M'C  Ifl  “  2  -  «  B  .  I 

1=1  It 

E!«  «l  !{„  ; 
eis  j  » ..  g 
“E?35^5ic  Zc  - 

fft  2  c  5  <  c  3  -o  .’5 

£.»  8  “ 

3=ui'“2?'B2#  ■ 

S53-2ilJ-8?“ 


1  ^  Ot~>  * 

JS'O  * 

-9  t?  « 

^  n  -  -  0 

-Wo—  c 

5  ^  >"+ 

1$  «*» 

^  OS  I  T 

XO-ci^ 

i  —  >r 
..  «■-  a  ; 


«  i  o-j  F  0*1 

ii-il*!' 

2  ^  0.  £.  q  ^  0 

[*  'T  E  ^ 

f-  o  1  f  2  u 

•  .:jS«’ 

* .;  r;  ^  = 

•  ■--•>:  E  O  > 

2  ^  i5  ,  U  O  O 

3?-5-s|^ 

^  ^  £  *  - 


■  Si  "  ; 

*5  Ex  ^  »*  < 

IT!^: 

c  J,  <9  ftj  «  c 

**  t:  ^  ^  -  « 
«  ?  <  C  3  *3 
O  j  .  «  0  O  1 
t.  v  ^  '*  • 

Sti£2-«  < 


5 

i2_;2  3 

-M3—  S 

»Sc’C  » 

*553,. 

if 

^  O  «  K  ^ 

!P 

,  a  v  H  ~ 

:i--g  .V 

1  .  £  "  -o 

>  <  g  V.  ST 


io;4  = 
i  m  ?  •»  ,0, 
-HO  2ii 


o  „  ri  .  ,  «  3 
Sc;S  -  «  c 

itow  w>  i>  •  P 
u  U  v  —  ..  0.  to  c  « 

uocuc^orjvij: 

'JirtOWSHHO* 

f  .  ■-  >.0  !r  g  c., 

.**  —  —  TJ  "  c  V  Q  2  t 

s«”  «o-2  £>“^v 

siiiNh  i%i 

T  S'  *j  ;**  S  «  ^  *  to  2 
x£^2*ofl,i5 

5  £h  :g^s  ui  o 

E  >■  “  2  ^  "  f 

s g »i~  ’■Sals 

=|I?!?Sg^S 

o.£  i*  £  f  2  ••* 

»e“«flS«eJS 

ssH^S^1 

S  ^  5  c  e,  t.  c  £  x 

JcJJSTo  _i! 

25  c  5r  "  «  r?  .S  T 

«  e  c  at  i 

u£*  3s  :J3^?S 

*=57(51>-'‘vc«a 

3  a  •-  I.  -.  3J  C 


